Low complexity decoding of low density parity check codes

ABSTRACT

An improved decoder and decoding method for low density parity check (LDPC) codes is provided. Decoding proceeds by repetitive message passing from a set of variable nodes to a set of check nodes, and from the check nodes back to the variable nodes. The variable node output messages include a “best guess” as to the relevant bit value, along with a weight giving the confidence in the guess. The check node output messages have magnitudes selected from a predetermined set including neutral, weak, medium and strong magnitudes. The check node output messages tend to reinforce the status quo of the input variable nodes if the check node parity check is satisfied, and tend to flip bits in the input variable nodes if the check node parity check is not satisfied. The variable node message weights are used to determine the check node message magnitudes.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 11/064,767, filed on Feb. 24, 2005 now U.S. Pat. No. 7,441,178 and entitled “Low Complexity Decoding of Low Density Parity Check Codes”.

FIELD OF THE INVENTION

This invention relates to decoding of low density parity check codes.

BACKGROUND

Coding is often used to reduce communication errors by deliberately introducing redundancy into a transmitted signal. When the signal is received, the redundancy introduced by the code can be used to detect and/or reduce errors. For example, a simple parity check code is obtained by transmitting blocks of N+1 bits, where N bits are data bits and one bit is a parity bit selected to make the parity of each N+1 bit block even. Such a code can provide detection, but not correction, of single bit errors. Introduction of more than one parity bit can improve code error reduction performance (e.g. by providing detection and/or correction of multiple bit errors). This code is an example of a block parity check code.

Block parity check codes can be considered more systematically in terms of a parity check matrix H. The matrix H has R rows and C columns, where C>R. Transmitted code words x are in the null space of H (i.e., Hx=0). Thus the columns of H correspond to symbols in the code word x (typically binary bits), and each row of H corresponds to a parity check condition on the code word x. Since a transmitted code word has C bits subject to R linear conditions, the data content of a code word is C−R bits if the rows of H are linearly independent. In some cases, the rows of H are not linearly independent, and in these cases the data content of a block is C−R*, where R* is the number of linearly independent rows of H (i.e., the dimension of the row space of H). When the rows of H are not linearly independent, H is transformed to an equivalent matrix H_(enc) having linearly independent rows for encoding. However, the original H matrix is still used for decoding. The rate of a block code is the ratio (C−R*)/C, and is a measure of the amount of redundancy introduced by the code. For example, a rate ½ code has one parity bit for each data bit in a block, and a rate ¾ code has one parity bit for each three data bits in a block.

A parity check code is completely defined by its parity check matrix H. Accordingly, encoding can be regarded as the process of mapping a sequence of data bits to code words in the null space of H. This encoding is typically done by constructing a generator matrix G from H such that a message vector u is mapped into a code word x in the null space of H via x^(T)=u^(T)G. Methods for constructing G given H are known in the art. For example, if H has the form [A|I] where A has dimensions n−k by k and I is an n−k dimensional identity matrix, G has the form [I|−A]. If H does not have this special form, G can still be constructed, but will not have the form [I|−A]. Similarly, decoding can be regarded as the process of estimating which code word was transmitted, given a received code word x′ which need not be in the null space of H due to transmission errors. Various methods for efficiently performing these encoding and decoding operations in practice have been developed over time.

In the course of this development, low density parity check (LDPC) codes have emerged as an especially interesting class of codes. The defining characteristic of an LDPC code is that the parity check matrix H is sparse (i.e., is mostly zeros). It is customary to use the notation LDPC(B, D) to refer to an LDPC code, where B is the total number of bits in a block, and D is the number of data bits in a block. Thus such a code has a parity check matrix H having B columns and B−D rows, if the rows are linearly independent. Some LDPC codes are referred to as “regular” codes because they have the same number d_(c) of non-zero elements in every row of H and have the same number d_(v) of non-zero elements in every column of H. Such codes are often referred to as (d_(v), d_(c)) LDPC codes. For example, a (3, 6) LDPC code has d_(v)=3 and d_(c)=6. In some cases, further structure has been imposed on H in order to improve encoding and/or decoding efficiency. For example, it is generally preferred for no two rows (or columns) of the H matrix to have more than one “1” in common.

The structure of regular LDPC codes can be appreciated more clearly in connection with a graph, as shown on FIG. 1. In the representation of FIG. 1, a set of variable nodes 110 and a set of check nodes 120 are defined. Each variable node is connected to d_(v) check nodes, and each check node is connected to d_(c) variable nodes. In the example of FIG. 1, d_(v)=3, d_(c)=6, and the connections from variable nodes to check nodes are not completely shown to preserve clarity. There is one variable node for each bit in a code word (i.e., there are C variable nodes), and there is one check node for each parity check condition defined by H (i.e., there are R check nodes). It is useful to define N(m) as the set of variable nodes connected to check node m, and M(n) as the set of check nodes connected to variable node n.

LDPC decoding can be regarded as a process of estimating values for the variable nodes given received variable data (which may have errors) subject to parity check conditions defined by each check node. Two approaches to decoding have been extensively considered: hard decision decoding and soft decision decoding. In hard decision decoding, received variable data is quantized to either of two binary values, and then error checking defined by the check nodes is performed on the quantized values. In the context of LDPC decoding, this approach includes “bit-flipping” decoding methods and majority-logic decoding methods. Suppose that a received code word has only a single bit error at variable node k. In this case, the check node conditions will be satisfied at all check nodes except for check nodes M(k). Since variable node k is the common element to all check nodes showing a violation, it should be flipped. While variations of such bit-flipping methods have been developed, a common feature of these approaches is quantization of received bit data to binary levels, followed by error correction processing.

As might be expected, such quantization of received bit data incurs a performance penalty, because information is lost. For example, if an analog signal between 0 and 1 is quantized to binary values of 0 and 1, received values of 0.51 and 0.99 are both quantized to 1. Clearly the “1” resulting from quantization of 0.51 is significantly less certain than the “1” resulting from quantization of 0.99. This performance penalty can be avoided by soft decision decoding. For LDPC codes, soft decision decoding is typically implemented as a message passing belief propagation (BP) algorithm. In such algorithms, variable messages are calculated in the variable nodes and passed to the check nodes. Next, check messages are computed in the check nodes and passed to the variable nodes. In both steps, computation of outgoing messages depends on received message inputs. These steps are repeated until a convergence condition is met (or a maximum number of iterations is reached). Soft decision decoding with a BP algorithm typically provides good performance, but such approaches tend to be more resource-intensive than hard decision approaches. Equivalently, soft decision decoding typically cannot be performed as rapidly as hard decision decoding.

This trade-off has motivated development of hybrid approaches having aspects of both hard and soft decision decoding to provide improved performance with reduced complexity. For example, weighted bit flipping (WBF) and modified WBF methods are described by Zhang and Fossorier in IEEE Comm. Lett., v8 n3, pp. 165-167, 2004. In these methods, bits are flipped responsive to calculations performed on unquantized received bits. These bit-flipping methods (as well as simple bit-flipping) typically require a search to be performed over all variable nodes. For example, in simple bit-flipping let ERR(n) be the number of parity check violations in the set of check nodes M(n) associated with variable node n. Bits corresponding to variable nodes having a maximal (or above-threshold) value of ERR(n) are flipped. In this example, a search is required to identify these variable nodes. Such a search can undesirably increase the computational resources required for decoding.

Hybrid decoding of LDPC codes is also considered in US patent application 2004/0148561, where various hybrid methods including both a bit flipping step and a belief propagation step are considered. Approaches of this kind do not provide simple implementation, since both bit-flipping and belief propagation are performed. In particular, these approaches do not alleviate the above-noted disadvantage of reduced decoding speed exhibited by conventional BP decoding.

Accordingly, it would be an advance in the art to provide hybrid decoding of LDPC codes that overcomes these disadvantages. More specifically, it would be an advance to provide a hybrid decoding method that does not require searching nodes for maximal (or above threshold) values and does not require implementation of conventional belief propagation.

SUMMARY

The present invention provides an improved decoder and decoding method for low density parity check (LDPC) codes. Decoding proceeds by repetitive message passing from a set of variable nodes to a set of check nodes, and from the check nodes back to the variable nodes. The variable node output messages include a “best guess” as to the relevant bit value, along with a weight giving the confidence in the guess (e.g., weak, medium or strong). The check node output messages have magnitudes selected from a predetermined set including neutral, weak, medium and strong magnitudes. The check node output messages tend to reinforce the status quo of the input variable nodes if the check node parity check is satisfied, and tend to flip bits in the input variable nodes if the check node parity check is not satisfied. The variable node message weights are used to determine the check node message magnitudes.

This method can be regarded as a hybrid method between conventional bit-flipping methods and conventional belief propagation (BP). More specifically, the underlying logic of the method of the present invention is bit flipping combined with the use of information (i.e., the variable message weights) giving the confidence in the bits being processed. The present invention may appear similar to BP because a message passing algorithm is employed, but the messages in the present invention are significantly simpler than in conventional BP. Furthermore, the present invention requires much less computation in the nodes (especially the check nodes) than in conventional BP.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows variable nodes and check nodes of an LDPC code as known in the art.

FIG. 2 shows a method for decoding an LDPC according to an embodiment of the invention.

FIG. 3 shows examples of LDPC(1024, 833) decoding performance, including an embodiment of the invention, for BPSK modulation.

FIG. 4 shows examples of LDPC(2048, 1723) decoding performance, including an embodiment of the invention, for BPSK modulation.

DETAILED DESCRIPTION

FIG. 2 shows a method for decoding LDPC codes according to an embodiment of the invention. This method can be appreciated as an improved bit-flipping algorithm, where quantized messages are passed to carry more information than in conventional bit-flipping methods. In addition, the quantization of the messages provides significant advantages of implementation simplicity compared to conventional belief propagation decoding by message passing. These aspects of the method of FIG. 2 will be further emphasized in the following description.

Step 202 on FIG. 2 is defining a set of variable nodes and a set of check nodes (e.g., as shown on FIG. 1). As indicated above, there is one variable node for each symbol (typically a bit) in a code block, and there is one check node for each parity check condition (i.e., for each row of the parity check matrix H). Since code performance tends to improve as block size increases, blocks are typically quite large (e.g., several thousand bits or more). In practice, the advantages of long blocks are traded off against disadvantages such as increased decoder latency to select a block size. The variable nodes and check nodes are indexed by integers n and m respectively. The function M(n) returns the set of check nodes associated with variable node n (i.e., the parity checks relating to variable node n), and the function N(m) returns the set of variable nodes associated with check node m (i.e., included in the parity check condition for check node m). LDPC codes suitable for decoding according to the invention can be regular or irregular. Regular codes have the same number of variable nodes associated with each check node, and have the same number of check nodes associated with each variable node. Equivalently, regular codes have the same number of ones in every row of H and the same number of ones in every column of H. Irregular codes have a variable number of ones in the rows and/or columns of H. The decoding methods of the present invention are generally applicable to LDPC codes (i.e., to any code having a sparse parity check matrix H). Thus H can be any random sparse matrix. The invention is also suitable for decoding LDPC codes having special properties (i.e., additional mathematical constraints on H). Such constraints can be useful in various ways, such as facilitating encoding.

Step 204 on FIG. 2 is predetermining neutral, weak, medium and strong check message magnitudes CK^(n), CK^(w), CK^(m), and CK^(s) respectively. According to the invention, the magnitude of the check messages is quantized to be one of these predetermined magnitudes, as discussed in more detail below. Such quantization of the check messages differs significantly from conventional belief propagation, where message quantization may occur incidentally (e.g., in a digital implementation) but is not a fundamental feature of the method.

Step 206 on FIG. 2 is initializing the variable nodes. Let W_(n) be the value of variable node n. At any point in the decoding process, the value W_(n) is an estimate of the received symbol corresponding to variable node n. It is preferred to represent these values as log-likelihood ratios (e.g., W_(n)=log(P₀(n)/P₁(n)), where P₁(n) is the probability that bit n is a one and P₀(n) is the probability that bit n is a zero. The log can be taken to any base, but the natural logarithm (i.e., base e) is preferred. The log-likelihood ratio can also be defined in terms of P₁(n)/P₀(n), naturally. The values W_(n) can also be expressed in terms of probabilities (e.g., W_(n)=P₀(n) or W_(n)=P₁(n)), but this representation is less preferred since it tends to require more computation to decode than the log-likelihood representation. Thus, the values W_(n) are preferably initialized to log-likelihood ratios L_(n) ⁽⁰⁾ calculated from received symbols. For example, large positive values of L_(n) ⁽⁰⁾ correspond to a high probability of a received 1, large negative values of L_(n) ⁽⁰⁾ correspond to a high probability of a received 0, and L_(n) ⁽⁰⁾=0 if the received symbol value is at the decision threshold between 0 and 1. As another example, for transmission of −1 or +1 over a channel having additive Gaussian noise having zero mean and variance σ², the log-likelihood ratio L_(n) ⁽⁰⁾ is 2y/σ² where y is the received symbol value. Suitable methods for computing L_(n) ⁽⁰⁾ from received symbol values are well known in the art.

Optionally, the initial values L_(n) ⁽⁰⁾ can be scaled by a scaling factor prior to performing the remaining steps of the decoding method. Such scaling can be regarded as a way to adjust decoding performance to match particular characteristics of the channel and/or modulation scheme being employed. As such, it is within the skill of an art worker to routinely adjust such scaling (e.g., by numerical simulation) to improve decoding performance. A scaling factor which is too large will tend to lead to overflow with associated loss of information. A scaling factor which is too small will tend to lose information by underflow or reduced precision. In typical receivers, automatic gain control is implemented to ensure input signals are within a well-defined range. Establishment of such a well-defined input signal range is helpful in selecting a suitable scaling factor for the initial values L_(n) ⁽⁰⁾. In practicing the invention, suitable scaling factors tend to be somewhat larger than in conventional decoding approaches. This difference is attributed to the weighting of variable messages as weak, medium, and strong (as discussed in detail below), and the high relative importance of weak variable messages.

Step 208 on FIG. 2 is quantizing the variable node values W_(n). More specifically, quantization is performed relative to a decision threshold X such that a quantized value Q_(n) is set to a predetermined value Q⁰ if W_(n)<X and is set to a predetermined value Q¹ if W_(n)≧X. For example, in the preferred case where the values W_(n) are log-likelihood ratios, a suitable decision threshold X is 0, and suitable values for Q⁰ and Q¹ are −1 and 1 respectively. Thus the quantized value Q_(n) can be regarded as the result of a hard decision. However, as indicated above, information is lost if a quantity such as W_(n) (which can take on a range of values) is replaced with a single bit of information (e.g., Q_(n)).

According to the invention, this issue is addressed by also computing a weight R_(n) to associate with the value Q_(n). The weight R_(n) takes on different values depending on how close the value W_(n) is to the threshold X. If W_(n) is close to X, the weight R_(n) is “weak”, since there is a relatively high probability that the corresponding value Q_(n) is incorrect. If W_(n) is far away from X, the weight R_(n) is “strong”, since there is a relatively low probability that the corresponding value Q_(n) is incorrect. For intermediate cases, the weight R_(n) is “medium”. More precisely, two weight thresholds T₁ and T₂ are defined, and R_(n) is weak if |W_(n)−X|≦T₁, R_(n) is medium if T₁<|W_(n)−X|≦T₂, and R_(n) is strong if T₂<|W_(n)−X|.

Step 210 on FIG. 2 is passing messages Z_(n) from the variable nodes to the check nodes. Each variable node n passes the same message Z_(n) to each of its associated check nodes M(n). The variable messages Z_(n) include the quantized values Q_(n) and the weights R_(n). Any mathematically unambiguous representation of Q_(n) and R_(n) can be used to construct the messages Z_(n). For example, weak, medium and strong weights can correspond to R_(n) equal to 1, 2, and 3 respectively, and Q_(n) can take on values of −1 or +1 as above. In this representation, the messages Z_(n) can be the product of Q_(n) and R_(n), since no ambiguity occurs.

Step 212 on FIG. 2 is calculating check node output messages L_(mn) from each check node m to its associated variable nodes N(m). Broadly speaking, each check node message provides information to its variable nodes as to whether or not the variable node values should be changed (e.g., by tending to flip a bit) based on the parity check condition at the check node. More specifically, the check node output messages have magnitudes which are selected from the predetermined magnitudes CK^(n), CK^(w), CK^(m), and CK^(s). Strong check node messages are sent in cases where it is clear whether or not to flip a bit, weak check node messages are sent in cases where it is less clear whether or not to flip a bit, medium check node messages are sent in intermediate cases, and neutral check node messages are sent in cases where it is not clear whether or not a bit should be flipped. The determination of how clear it is whether or not to flip a bit depends on the weights R_(n) (i.e., confidence) of the received variable node messages Z_(n). Mathematically, the messages L_(mn) are determined by a parity check function f_(mn)[{Z_(n):nεN(m)}] of the variable messages provided to check node m.

In a preferred embodiment of the invention, the parity check function is computed according to the following method. In this preferred embodiment, we assume Q⁰=−Q¹ (e.g., Q⁰ can be 1 and Q¹ can be −1, or vice versa). First, a parity check condition for check node m is checked. More specifically, a determination is made as to whether or not the quantized values Q_(n) input to check node m satisfy the parity check condition for check node m. For example, a commonly employed parity condition is satisfied if and only if an even number of associated variable nodes have quantized values Q_(n) equal to one. Thus, if the parity check condition is satisfied, there are either no bit errors, or an even number (e.g., 2, 4, etc.) of bit errors. If the parity check condition is not satisfied, there is an odd number (e.g., 1, 3, etc.) of bit errors. In most cases, systems are designed to have a low probability of error, which means that if parity check is satisfied, it is most probable that all bits are correct, and if parity check is violated, it is most probable that only one bit is in error. These considerations will be apparent in the following steps.

Second, the check node confidence is established. More specifically, this step entails setting a weak count WK_(m) and a medium count MD_(m) equal to a count of the number of weak and medium weights respectively included in the input variable messages {Z_(n):nεN(m)} to check node m. If there is only one weak variable message, an index n_(weak) is set equal to the index n of the variable node in N(m) having a weak weight. Similarly, if there is only one medium variable message, an index n_(med) is set equal to the index n of the variable node in N(m) having a medium weight.

Finally, the parity check function f_(mn) is computed. Two cases are distinguished. If the parity check condition is satisfied, f_(mn) is computed according to

$\begin{matrix} {f_{mn} = \left\{ \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {\underset{\_}{\begin{matrix} {Q_{n}{CK}^{s}} & {{if}\mspace{14mu}\begin{matrix} {{{WK}_{m} = 0};} \\ {{{MD}_{m} = 0};} \\ {{{for}\mspace{14mu} n} \in {N(m)}} \end{matrix}} \end{matrix}}\mspace{259mu}} \\ {\underset{\_}{\left\{ {\begin{matrix} {{Q_{n}{CK}^{s}\mspace{14mu}{for}\mspace{14mu} n} = n_{med}} \\ {{Q_{n}{CK}^{m}\mspace{14mu}{for}\mspace{14mu} n} \neq n_{med}} \end{matrix}\mspace{14mu}{if}\mspace{14mu}\begin{matrix} {{{WK}_{m} = 0};} \\ {{{MD}_{m} = 1};} \\ {{{for}\mspace{14mu} n} \in {N(m)}} \end{matrix}} \right.}\mspace{121mu}} \end{matrix} \\ {\underset{\_}{{{Q_{n}{CK}^{w}\mspace{14mu}{if}\mspace{14mu}{WK}_{m}} = 0};{{MD}_{m} \geq 2};{{{for}\mspace{14mu} n} \in {N(m)}}}\mspace{65mu}} \end{matrix} \\ {\underset{\_}{\left\{ {\begin{matrix} {{Q_{n}{CK}^{s}\mspace{14mu}{for}\mspace{14mu} n} = n_{weak}} \\ {{Q_{n}{CK}^{w}\mspace{14mu}{for}\mspace{14mu} n} \neq n_{weak}} \end{matrix}\mspace{14mu}{if}\mspace{14mu}\begin{matrix} \begin{matrix} {{{WK}_{m} = 1};} \\ {{{MD}_{m} = 0};} \end{matrix} \\ {{{for}\mspace{14mu} n} \in {N(m)}} \end{matrix}} \right.}\mspace{121mu}} \end{matrix} \\ {\underset{\_}{{{Q_{n}{CK}^{w}\mspace{14mu}{if}\mspace{14mu}{WK}_{m}} = 1};{{MD}_{m} \geq 1};{{{for}\mspace{14mu} n} \in {N(m)}}}\mspace{40mu}} \end{matrix} \\ {{\begin{matrix} {CK}^{n} & {{{{if}\mspace{14mu}{WK}_{m}} \geq 2};{{{for}\mspace{14mu} n} \in {N(m)}}} \end{matrix}.}\mspace{124mu}} \end{matrix} \right.} & (1) \end{matrix}$ If the parity check condition is not satisfied, f_(mn) is computed according to

$\begin{matrix} {f_{mn} = \left\{ \begin{matrix} \begin{matrix} \begin{matrix} {\underset{\_}{\begin{matrix} {{- Q_{n}}{CK}^{w}} & {{{{if}\mspace{14mu}{WK}_{m}} = 0};{{{for}\mspace{14mu} n} \in {N(m)}}} \end{matrix}}\mspace{110mu}} \\ {\underset{\_}{\left\{ {\begin{matrix} {{{- Q_{n}}{CK}^{s}\mspace{14mu}{for}\mspace{14mu} n} = n_{weak}} \\ {{{- Q_{n}}{CK}^{w}\mspace{14mu}{for}\mspace{14mu} n} \neq n_{weak}} \end{matrix}\mspace{14mu}{if}\mspace{11mu}\begin{matrix} \begin{matrix} {{{WK}_{m} = 1};} \\ {{{MD}_{m} = 0};} \end{matrix} \\ {{{for}\mspace{11mu} n} \in {N(m)}} \end{matrix}} \right.}\mspace{76mu}} \end{matrix} \\ \underset{\_}{\begin{matrix} {{- Q_{n}}{CK}^{w}} & {{{{if}\mspace{14mu}{WK}_{m}} = 1};{{MD}_{m} \geq 1};{{{for}\mspace{14mu} n} \in {N(m)}}} \end{matrix}} \end{matrix} \\ \begin{matrix} {CK}^{n} & {{{{{if}\mspace{14mu}{WK}_{m}} \geq 2};{{{for}\mspace{14mu} n} \in {{N(m)}.}}}\mspace{59mu}} \end{matrix} \end{matrix} \right.} & (2) \end{matrix}$ Alternative forms for these equations are possible. For example, if the parity check condition is not satisfied, f_(mn) can also be computed according to

$\begin{matrix} {f_{mn} = \left\{ \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {\underset{\_}{\begin{matrix} {{- Q_{n}}{CK}^{w}} & {{if}\mspace{14mu}\begin{matrix} {{{WK}_{m} = 0};} \\ {{{MD}_{m} = 0};} \\ {{{for}\mspace{14mu} n} \in {N(m)}} \end{matrix}} \end{matrix}}} \\ {\underset{\_}{\left\{ {\begin{matrix} {{{- Q_{n}}{CK}^{s}\mspace{14mu}{for}\mspace{14mu} n} = n_{med}} \\ {{{- Q_{n}}{CK}^{w}\mspace{14mu}{for}\mspace{14mu} n} \neq n_{med}} \end{matrix}\mspace{14mu}{if}\mspace{14mu}\begin{matrix} {{{WK}_{m} = 0};} \\ {{{MD}_{m} = 1};} \\ {{{for}\mspace{14mu} n} \in {N(m)}} \end{matrix}} \right.}\mspace{121mu}} \end{matrix} \\ {\underset{\_}{{{{- Q_{n}}{CK}^{w}\mspace{14mu}{if}\mspace{14mu}{WK}_{m}} = 0};{{MD}_{m} \geq 2};{{{for}\mspace{14mu} n} \in {N(m)}}}\mspace{65mu}} \end{matrix} \\ {\underset{\_}{\left\{ {\begin{matrix} {{{- Q_{n}}{CK}^{s}\mspace{14mu}{for}\mspace{14mu} n} = n_{weak}} \\ {{{- Q_{n}}{CK}^{w}\mspace{14mu}{for}\mspace{14mu} n} \neq n_{weak}} \end{matrix}\mspace{14mu}{if}\mspace{14mu}\begin{matrix} \begin{matrix} {{{WK}_{m} = 1};} \\ {{{MD}_{m} = 0};} \end{matrix} \\ {{{for}\mspace{14mu} n} \in {N(m)}} \end{matrix}} \right.}\mspace{121mu}} \end{matrix} \\ {\underset{\_}{{{{{- Q_{n}}{CK}^{w}\mspace{14mu}{if}\mspace{14mu}{WK}_{m}} = 1};{{MD}_{m} \geq 1};{{{for}\mspace{14mu} n} \in {N(m)}}}\mspace{14mu}}\mspace{40mu}} \end{matrix} \\ {{\begin{matrix} {CK}^{n} & {{{{if}\mspace{14mu}{WK}_{m}} \geq 2};{{{for}\mspace{14mu} n} \in {N(m)}}} \end{matrix}.}\mspace{124mu}} \end{matrix} \right.} & \left( {2a} \right) \end{matrix}$ Some comparisons have been performed between decoding according to Eqs. 1 and 2 and decoding according to Eqs. 1 and 2a. In these comparisons, no significant performance difference was seen. It is helpful to regard the check node output messages as suggestions about how to modify the variable node values W_(n) based on the messages input to the check nodes. Consideration of specific cases relating to equations 1 and 2 will clarify these equations, and also suggest further alternatives and variations of such equations to the skilled art worker.

For example, if there are two or more weak variable messages at a check node, there is not enough confidence to suggest anything definite, so neutral messages CK^(n) are sent back. At the opposite extreme, if there are no weak or medium input messages and parity check is satisfied, everything appears to be in order and strong output messages Q_(n)CK^(s) are sent back reinforcing the status quo. If there are no weak input messages and parity check is violated, something is wrong, but there is no clear single target for bit flipping, so a weak bit flipping output message −Q_(n)CK^(w) is sent out. Note that a check message proportional to Q_(n) will act to reinforce the status quo and a check message proportional to −Q_(n) will act to change the status quo (by flipping a bit).

If parity check is satisfied, and there are no weak input messages and two or more medium input messages, everything is apparently in order, but with relatively low confidence. Therefore, weak output messages Q_(n)CK^(w) are sent out reinforcing the status quo. Such messages Q_(n)CK^(w) are also sent out if parity check is satisfied, there is one weak input message and one or more medium input messages. If parity check is violated and there is one weak input message and one or more medium input messages, weak bit-flipping output −Q_(n)CK^(w) messages are sent out.

If parity check is violated, and there is one weak input message and no medium-input messages, there is a clear indication of which bit should most probably be flipped. Accordingly, a strong bit-flipping message −Q_(n)CK^(s) is sent to the weak variable node (i.e., n=n_(weak)), and weak bit-flipping messages −Q_(n)CK^(w) are sent to the other variable nodes (i.e., n≠n_(weak)). Similarly, if parity check is satisfied and there is one weak input message and no medium input messages, a strong status quo message Q_(n)CK^(s) is sent to the weak variable node (i.e., n=n_(weak)) and weak status quo messages Q_(n)CK^(w) are sent to the other variable nodes (i.e., n≠n_(weak)).

Finally, if parity check is satisfied and there are no weak input messages and a single medium input message, a strong status quo message Q_(n)CK^(s) is sent to the medium variable node (i.e., n=n_(med)) and medium status quo messages Q_(n)CK^(m) are sent to the other variable nodes (i.e., n≠n_(med)).

Step 214 on FIG. 2 is updating the variable node values W_(n) responsive to the check messages L_(mn). In the preferred representation where the variable values W_(n) are log-likelihood ratios, this updating is calculated according to

$\begin{matrix} {W_{n} = {L_{n}^{(0)} + {\sum\limits_{m \in {M{(n)}}}^{\;}\;{L_{mn}.}}}} & (3) \end{matrix}$ Thus W_(n) is updated to be the sum of the initial estimate L_(n) ⁽⁰⁾ and all of the messages L_(mn) from the associated check nodes M(n). Here the effect of the check message magnitudes can be most clearly appreciated. For example, one possible assignment of check message magnitudes is CK^(n)=0, CK^(w)=1, CK^(m)=2 and CK^(s)=3. Thus a neutral message has no effect on a variable node it is sent to, and the effect of the other messages is larger for strong than for medium, and larger for medium than for weak.

In an alternative embodiment of the invention, the updating of variable node values is calculated according to

$\begin{matrix} {W_{n} = {{f\left( {P\; C\; C} \right)} + L_{n}^{(0)} + {\sum\limits_{m \in {M{(n)}}}^{\;}\; L_{mn}}}} & \left( {3a} \right) \end{matrix}$ where PCC is a parity check count and f(PCC) is a function of this parity check count. More specifically, PCC is the number of parity errors present in the set of check nodes M(n) associated with variable node n. For example, |f(PCC)| can be one if PCC is one or less, and f(PCC) can be zero if PCC is two or greater. The sign of f(PCC) is selected to reinforce the value of W_(n) (i.e., is the same as the sign of

$\left. {L_{n}^{(0)} + {\sum\limits_{m \in {M{(n)}}}^{\;}\; L_{mn}}} \right).$

In this manner, variable nodes having a small number of parity errors in their associated check nodes have their values increased, which tends to increase the confidence weights assigned to these variable nodes. The use of this f(PCC) in cases otherwise similar to the examples of FIGS. 3-5 has improved performance by about 0.15 dB. Other functions f(PCC) having the same general characteristics may be more suitable than this exemplary f(PCC) for certain cases, and it is within the skill of an art worker to determine such functions (e.g., by performing routine numerical simulations).

Steps 208, 210, 212 and 214 are repeated in sequence until a termination condition is satisfied. The termination condition can be satisfaction of the parity check condition at all check nodes. The termination condition can also be completion of a predetermined number of iterations of this sequence. A preferred termination condition is to end decoding if all parity check conditions are satisfied, or if a maximum number of iterations is reached, whichever happens first.

Several advantages of the present invention are evident from the preceding exemplary description of a preferred embodiment. In particular, no search is required of either the variable nodes or check nodes, which is a significant advantage compared to some conventional bit flipping methods. Furthermore, the decoding computations of the present invention are much simpler than the computations of conventional belief propagation decoding (especially in the check nodes).

FIGS. 3-4 show simulated decoding performance for various LDPC decoding methods, including embodiments of the present invention, for an additive white Gaussian noise channel. On these plots, SBF refers to the simple bit-flipping method, MWBF refers to modified weighted bit flipping as discussed above, and XWBF refers to decoding according to embodiments of the present invention. For these examples, decoding is performed according to Equations 1, 2, and 3. In these embodiments, the log-likelihood representation for W_(n) is used, and the decoding parameters are set as follows: CK^(n)=0, CK^(w)=1, CK^(m)=2, CK^(s)=3, Z=0, T₁=2, T₂=7, Q⁰=1, and Q¹=−1. Although these parameter values have been found suitable for these examples, adjustment of these values is likely to be necessary to obtain optimal results for various other channels and/or modulation schemes. Such adjustment is routine, and well within the skill of an art worker. For example, numerical simulations are suitable for optimizing numerical decoding parameters.

The weak message threshold T₁ is noteworthy. If it is too high, then too many variable messages will be classified as weak. In particular, every check node may have two or more weak input messages. Such a state is undesirable, since all check messages will be neutral, and decoding will therefore make no progress. Thus the threshold T₁ should be set low enough so that this does not occur. Alternatively, an input scaling factor can be increased, such that the same threshold T₁ is relatively smaller compared to typical variable node values.

Modulation formats having more than one bit of information per symbol can be used with binary LDPC codes in known ways. For example, a modulation format having four levels 0, 1, 2, and 3 (i.e., PAM-4) can be mapped into bit patterns 00, 01, 11, and 10 respectively, which can be encoded, transmitted and decoded using a binary LDPC code. On FIGS. 3-4, the vertical axis is bit error rate, and the horizontal axis is E_(b)/N₀, where E_(b) is the bit energy and N₀ is the noise power spectral density. The examples of FIGS. 3-4 all relate to binary phase shift keying (BPSK) modulation.

FIG. 3 shows LDPC(1024,833) decoding performance for BPSK modulation, and FIG. 4 shows LDPC(2048, 1723) decoding performance for BPSK modulation. The codes for FIGS. 3 and 4 are regular LDPC codes. For the LDPC(1024,833) code of FIG. 3, dv=10 and dc=32. For the LDPC(2048,1723) code of FIG. 4, dv=6 and dc=32. All decoding methods were computed with 3-bit fixed point quantization of log-likelihood ratios. In both examples, the performance of the present invention (i.e., XWBF) is better than the bit flipping methods (i.e., more than 1 dB better than SBF and more than 0.7 dB better than MWBF) Thus an advantage of the present invention is to provide decoding performance that is better than conventional bit flipping methods without being computationally demanding.

The preceding description is by way of example as opposed to limitation. Accordingly, many variations of the above examples also fall within the scope of the invention. For example, the preceding description is in terms of log-likelihood ratios. A mathematically equivalent formulation in terms of probabilities (as known in the art for conventional LDPC decoding) can also be employed to practice the invention.

In the preceding description, inequalities were used to define ranges. For completeness, the case of mathematical equality was explicitly provided for (e.g., as in “R_(n) is medium if T₁<|W_(n)−X|≦T₂”). The seemingly different “R_(n) is medium if T₁≦|W_(n)−X|<T₂” is essentially equivalent to the preceding quote. Thresholds such as T₁ and T₂ will typically be set by a process of numerical experimentation, which can be performed with any consistent arrangement of the equality conditions.

The preceding description has focused on decoding methods according to the invention. A decoder according to the invention includes a processor which carries out the steps of a method of the invention. Such a processor can be implemented in any combination of hardware and/or software. Numerical calculations in such a processor can be performed in a fixed point representation and/or in a floating point representation. A preferred embodiment of a decoder of the invention is implemented in VLSI circuitry to maximize decoding speed. 

1. A method for decoding a check code, the method performed by circuitry and comprising: quantizing a value stored at each of a plurality of variable nodes to generate a plurality of respective quantized values, wherein the value stored at each of the plurality of variable nodes is based, at least in part, on a received symbol; assigning a weight to each of the respective quantized values, wherein each weight is based, at least in part, on a difference between the stored value and the quantized value; transmitting a variable message from each of the plurality of variable nodes to a respective associated plurality of connected check nodes, wherein each variable message is based, at least in part, on the quantized value and weight associated with the respective variable node; calculating a check node output message at each of the plurality of check nodes, wherein each check node output message has a magnitude based, at least in part, on the weight associated with the quantized value received by the respective check node; and updating at least one stored value in response to the check node output messages.
 2. The method of claim 1, wherein each of the plurality of variable nodes corresponds to a received symbol.
 3. The method of claim 1, wherein each of the plurality of check nodes corresponds to a parity check condition.
 4. The method of claim 1, further comprising initializing the plurality of variable nodes prior to said quantizing a value stored at each of a plurality of variable nodes.
 5. The method of claim 4, wherein said initializing the plurality of variable nodes comprises initializing each variable node to equal a log-likelihood ratio calculated from a received symbol.
 6. The method of claim 5, wherein said initializing the plurality of variable nodes further comprises scaling the log-likelihood ratios by a scaling factor.
 7. The method of claim 1, wherein said assigning a weight to each of the respective quantized values comprises assigning a weight to each of the respective quantized values that is selected from a group consisting of a weak weight, a strong weight, and a medium weight.
 8. The method of claim 7, wherein said assigning a weight to each of the respective quantized values further comprises: assigning the weak weight if a difference between the respective quantized value and the value stored at the variable node is less than or equal to a first threshold; assigning the medium weight if a difference between the respective quantized value and the value stored at the variable node is greater than the first threshold but less than or equal to a second threshold; and assigning the strong weight if a difference between the respective quantized value and the value stored at the variable node is greater than the second threshold.
 9. The method of claim 1, wherein said transmitting a variable message from each of the plurality of variable nodes comprises transmitting a same message from at least one of the plurality of variable nodes to each of the plurality of check nodes associated with the at least one variable node.
 10. The method of claim 1, wherein each respective variable message comprises a product of the weight and the quantized value associated with the respective variable node.
 11. The method of claim 1, wherein said calculating a check node output message comprises selecting a magnitude of each check node output message from a group of predetermined check node message magnitudes consisting of a neutral magnitude, a weak magnitude, a medium magnitude, and a strong magnitude.
 12. The method of claim 1, wherein said calculating a check node output message comprises: checking a parity check condition for each of the plurality of check nodes; establishing a confidence for each of the plurality of check nodes, wherein the confidence for each check node is based, at least in part, on the weights of the variable messages received by the respective check node; and computing a parity check function for each of the check nodes based, at least in part, on the confidence of the respective check node.
 13. The method of claim 12, wherein the weights are selected from a group consisting of a weak weight, a medium weight, and a strong weight, and wherein said computing a parity check function is based, at least in part, on a number of variable messages received by the respective check nodes having a weak weight.
 14. The method of claim 12, wherein said computing a parity check function comprises computing the parity check function according to a first set of conditions if the parity check condition is satisfied and computing the parity check function according to a second set of conditions if the parity check condition is not satisfied.
 15. The method of claim 12, further comprising repeating transmitting a variable message to the respective check nodes, calculating check node output messages, and updating the variable node stored values until a termination condition is satisfied.
 16. A decoder including circuitry comprising: a plurality of variable nodes, each variable node configured to store a respective value corresponding to a portion of a parity check code and generate a variable message, wherein the variable message includes a quantized value corresponding to the respective value and a weight based, at least in part, on a difference between the respective value and the quantized value; and a plurality of check nodes, each coupled to a respective set of the plurality of variable nodes and configured to receive the variable messages generated by the respective set of the plurality of variable nodes, wherein each of the plurality of check nodes is further configured to generate a check node output message having a magnitude based, at least in part, on the quantized values received by the respective check node and transmit the check node output message to the respective set of variable nodes, and wherein each of the respective set of variable nodes is configured to update the respective value in response to the check node output message.
 17. The decoder of claim 16, wherein the updated values stored at the variable nodes represent an estimate of transmitted symbols.
 18. The decoder of claim 16, wherein each of the plurality of variable nodes corresponds to a received symbol.
 19. The decoder of claim 16, wherein each of the plurality of check nodes corresponds to a parity check condition.
 20. The decoder of claim 16, wherein each value stored at the plurality of variable nodes is initially equal to a log-likelihood ratio calculated from a received symbol prior to being updated.
 21. The decoder of claim 20, wherein the log-likelihood ratios are scaled by a scaling factor.
 22. The decoder of claim 16, wherein the weight of each of the respective quantized values is selected from a group consisting of a weak weight, a medium weight, and a strong weight.
 23. The decoder of claim 22, wherein the weak weight is used if a difference between the respective quantized value and the value stored at the variable node is less than or equal to a first threshold, the medium weight is used if a difference between the respective quantized value and the value stored at the variable node is greater than the first threshold but less than or equal to a second threshold, and the strong weight is used if a difference between the respective quantized value and the value stored at the variable node is greater than the second threshold.
 24. The decoder of claim 16, wherein at least one of the variable nodes is further configured to transmit a same message to each of the plurality of check nodes associated with the at least one variable node.
 25. The decoder of claim 16, wherein the variable message generated by each variable node comprises a product of the weight and the quantized value associated with the respective variable node.
 26. The decoder of claim 16, wherein the check node output message is generated, at least in part, by selecting a magnitude of the check node output message from a group of predetermined check node message magnitudes consisting of a neutral magnitude, a weak magnitude, a medium magnitude, and a strong magnitude.
 27. The decoder of claim 16, wherein each of the plurality of check nodes are further configured to generate the check node output message, at least in part, by: checking a parity check condition for the check node; establishing a confidence for the check node, wherein the confidence is based, at least in part, on the weights of the variable messages received by the check node; and computing a parity check function for the check node based, at least in part, on the confidence of the check node.
 28. The decoder of claim 27, wherein the weights are selected from a group consisting of a weak weight, a medium weight, and a strong weight, and wherein the plurality of check nodes are further configured to compute the parity check function based, at least in part, on a number of variable messages received by the check node having a weak weight.
 29. The decoder of claim 27, wherein the plurality of check nodes are further configured to compute the parity check function according to a first set of conditions if the parity check condition is satisfied and to compute the parity check function according to a second set of conditions if the parity check condition is not satisfied.
 30. The decoder of claim 27, wherein the variable nodes are further configured to repeatedly transmit a variable message to the respective check nodes and update the variable node stored values until a termination condition is satisfied. 